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AT LAKEVIEW 

ABSTRACT 

The cooling water plumes from an operating 
thermal electric generating station at Lakeview, west 
of Toronto, on Lake Ontario, were defined by water 
temperature measurements. The generating station 
discharged between 1400 to 2500 cfs of cooling water 
at 15 °F above intake water temperature. The temperature 
measurements were made throughout the year at various 
depths in proximity of the discharge to define the 
physical cooling water plume in the lake. A method for 
estimating temperature isotherms was developed by 
statistically analyzing the temperature measurements. 
Equations were developed which defined pliame area as a 
function of distance from the discharge, temperature as 
a function of plume area and plume centreline as a 
function of wind. Different relationships were developed 
for summer and winter conditions based upon the ambient 
lake temperature. The most significant cooling resulted 
from mixing with lake water. 



ESTIMATE OF THE PHYSICAL 
EXTENT OF THERMAL PLUMES 
AT LAKEVIEW 

INTRODUCTION 

With the advent of large thermal electrical 
generating stations employing Great Lakes water for cooling, 
there has been increasing concern on the possible detrimental 
effects from the addition of heat to the water environment. 
Research on the effect of heat inputs on the water environ- 
ment has been extensive (Parker, 1969 and ASCE, 1967). 
However, the methods for predicting the physical extent of 
cooling water plumes on lakes are notably few. This aspect 
is considered essential by the Ontario Water Resources 
Commission. Once defined, the local effects of the heated 
plume on the water environment could be examined. 

Studies were undertaken by the Hydro-Electric Power 
Commission of Ontario (H.E.P.C.) to define the physical extent 
of the cooling water plume at the Lakeview generating station 
on Lake Ontario. Lakeview discharged approximately 1400-2500 
cfs at 15°F above intake water temperature. Initially, the 
studies consisted of measuring water temperatures with a 
telethermometer at predetermined locations in proximity of 
the discharge throughout the year. Normally, a survey took 
less than 3 hours. Meteorological data collected during the 
surveys consisted of wind and air temperature determinations . 
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This report discusses the development of a 
first order estimate of the spatial distribution of 
isotherms within the cooling water plume. The estimate 
is a first order approximation based on limited survey 
data. A better estimate could be made by solving the 
heat budget equation. However, the solution of the heat 
budget equation requires more extensive measurements of 
wind, relative humidity, air and water temperatures, 
net solar radiation and water movement patterns. 
Measurements of these parameters are now being made by 

H.E.P.C. 

Csanady (1970) recently showed that rate of 
heat lost to atmosphere is one-third of that lost to the 
lake in studies at Douglas Point on Lake Huron. Palmer 
(1969) also showed that excess heat in a plume is 
dissipated through horizontal mixing with the main lake. 
It will be shown that the principal cooling mechanism 
appears to be mixing with the ambient lake water. These 
results are in conflict with Brady (1969) and Parker (1969) 
who both indicate that the main cooling is atmospheric. 
Rieshol (1971) in studies for a cooling water discharge on 
Dardanelle Lake in Arkansas found that cooling by entrainment 
of lake water is the main mechanism in the first several 
thousand feet while atmospheric cooling prevails farther away 
from the discharge point. 
EXPERIMENTAL METHOD 

Water temperature surveys were conducted by H.E.P.C. 
with a small fast boat on relatively calm days. Surface 



and some depth temperatures were measured at the 
intersection points of a 1000 feet square geometric 
grid. The intersection points were defined by visual 
alignment on shore range poles with some assistance 
from shore based optical transits when required. All 
temperatures were measured with a Yellow Springs thermister 
telethermometer having an accuracy of 0.2 F. 

An average water temperature survey was 
completed in less than three hours. During the period 
of survey, the following parameters other than water 
temperature were measured several times on shore and 
then averaged. 

1. Air temperature in "F 

2 . Wind speed in mph 

3. Wind direction in degrees from north 

4. Rate of discharge from outfall in cfs 
Observed isotherms were plotted for each survey. The 
centreline of the plxime was drawn by joining the farthest 
points of each isotherm from the outfall. The area 
enclosed by each isotherm was measured by planimeter. 

PRELIMINARY INVESTIGATIONS 

Several different methods for predicting lake 
temperature isotherms from the available data were tried 
with varying degrees of success in terms of matching actual 
observed isotherms. Some of the methods tried are described 
briefly. 
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Gifford's Normal Distribution Dispersion Model 
This method assumes that the concentration 
distribution of any continuously discharged material in a 
moving stream will be normally distributed across any 
plane perpendicular to the moving stream. As the form of 
the concentration distribution is defined, i.e. a normal 
distribution, it is only necessary to define spread 
(variance) and current as functions of downstream location 
to determine concentrations at different distances from the 
source. This method has been applied successfully in 
meteorology (Gifford, 1955), ocean discharges (Foxworthy, 
1966) and lakeshore areas (Palmer, 1969). 

? (x) = 2Q (1) 

where 

C (x) = maximum mean concentration at distance 
max 

X from source in lbs/ft* 
Q = continuous mass flow in lbs/sec 
U = current velocity in ft/sec 
o^{o] = variance of concentration distribution 

at discharge in ft^ 
0^{x} = variance of concentration distribution 

at distance x from discharge point in ft 
This method was applied with some success to a cooling water 
plume in May 1968 (Palmev, 1969). The method considers 
dilution with ambient lake water. Since lake cooling has 
been shown as the principal cooling mechanism at least in the 
first few thousand feet, this method could be used effectively 



provided current velocities are known. Current velocities 
in near shore areas with shore topographical effects, could 
not be satisfactorily predicted from wind velocities. 
Therefore, the method was rejected. 

Johns Hopkin's Method 

This method is based on statistical fitting of the 
principle mechanisms responsible for cooling in small ponds 
or lakes (Brady, 1969). The technique was successfully 
applied to experimentally determined plumes with the computed 
parameters comparing favourably with other studies. However, 
it was necessary to define the portion of the lake affected 
by the plume to apply the method. For small lakes, this is 
simply the total surface area. As it was not possible to 
satisfactorily define the area of the lake affected, this 
method was abandoned. 

PROPOSED METHOD 

The above methods described briefly could not be 
used successfully with the available survey data. Therefore, a 
statistical approach was used to develop a prediction model 
for the thermal plume. The data used is listed in Appendix 2. 

The field survey data was divided into two seasons, 
i.e. winter and summer since the entire data set produced 
statistically weaker equations with large standard errors of 
estimates. The equations were improved by separating the data 
into 'winter' and 'summer' sets. Since much of cooling is due 
to mixing with lake water and air temperatures are very 



variable, ambient lake water temperature was chosen as 
the criterion for separation of data. Lake temperatures 
of 42°F, 45°F and 47''F were tried for separation of data. 
45 °F was found to be the best since it produced 
statistically best prediction equations for both sets. 
Best fit equations were developed for the following: 

a. area (A) contained by isotherms and the 
distance (X) (the farthest point on isotherm 
from the outfall, measured along the plume 
centreline) 

b. area (A) and temperature difference (ATp^^) 
between isotherm corresponding to area and 
ambient lake temperature 

c. distance (X) and temperature difference (ATp^^) 
From the first two relations, it is possible to 

define plume area and temperature at any distance (X) from 
the outfall. These equations are necessary to determine 
the plume temperature and area at any distance from outfall 
to enable plotting of the plume. The relation between X 
and AT_,^ has been obtained for completeness, but cannot be 

PL 

used to determine isotherm shapes since it does not define 
isotherm area. Various shapes, namely circular, elliptical 
and one leaf of multi-leaved rose were investigated as 
possible shapes for the predicted isotherm. One leaf of 
multi-leaved rose was found to approach the observed 
isotherm shape best and was adopted in the present model. 
The direction of the plume centreline was determined frcxn 
the wind vector. This method will be described in the 
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section 'Definition of Plume Centreline'. 

Plume Area and Distance 

Relations have been found for area or spread or 
variance as a function of distance from a source by many 
experimenters lAlsafar^ 1966; Csanady, 1964; Okubo, 1968; 
Foxworthyj 1966 and Murthy , 1969). In the case of 
near-shore areas of the Great Lakes, the relations appear 
to be general in slope for a particular location, but vary 
in intercept for the maximum, mean and minimum conditions 
(Palmer and Izatt, 1970). On this basis, it is probable 
that a statistically significant equation could be 
determined relating area (contained by the isotherms) of 
the plume and distance from outfall along the plume centre- 
line. The best fitted equation (see Table 1 and Figure 
2) for the summer data set was : 

Area = -0.696 + 1.384 (Distance) (2) 

with F, Qc = 117.1 : multiple correlation 
coefficient - 0.74: standard error - 3.11: 
standard deviation area = 4.62 
The best winter data equation (see Figure 3) was: 

Area = -4.047 + 2.43 (Distance) (3) 

with F, ,, = 283.6: multiple correlation 

J. , o J 

coefficient = 0.90: standard error = 2.49: 

standard deviation area = 5.80. 
It is immediately apparent that there was a statistically 
significant relationship between the areas contained by the 
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isotherms and the distance from the outfall. The 
winter equation was statistically stronger, when 
lake temperature variations were smallest. 

Isotherm Area and Temperature Difference (Plume and Ambient 
Lake Temperature) 

The cooling of the plume to either the 
atmosphere or the lake must be related to the plume area. 
The entrainment of lake water at ambient temperature is 
in part a function of the plume area in contact with the 
lake water. Similarly, the amount of atmospheric cooling 
is in part a function of the surface area. If one assumes 
that the plume area in contact with the atmosphere is 
nearly the same as that in contact with the lake, one 
would expect a relationship to exist between the area 
contained by an isotherm and the temperature difference 
between plume and ambient lake temperature. This seems 
to be a valid assumption as most of the thermal pliime 
appears to be confined to the surface regions of the lake, 
A statistically significant relationship was found (see 
Table 3 and Figure 4) . The best equation for the summer 
data including a flow term was: 

Area = 90.53 - 8.86 In (flow X temperature 

difference) (4) 

with F, Qi- = 135.2: multiple correlation 
coefficient = 0.77: standard error = 2.99: 
standard deviation area = 4.62. 
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The best winter equation including a flow term {see Figure 5) was: 
Area = 74.0 - 6.88 In (flow X 
temperature difference) (5) 

with F, ^-, = 126.3: multiple correlation 

coefficient = 0.82: standard error = 3.37: 

standard deviation area = 5.80. 

The question immediately arises whether a better 
equation could be determined if temperature differences 
between the plume and atmosphere as well as the temperature 
difference between ambient lake and plume were included. 
Statistical fitting of the summer and winter data 
including both temperature differences is tabulated in 
Table 2, Significant equations were found relating the 
temperature differences and isotherm area and for the 
summer data the standard error was slightly less. 
However "t" tests on the significance of the two temper- 
ature differences indicate that the air-plume difference 
was not significant for the summer data and not significant 
at a= 0.01 for the winter data. The predominate cooling 
mechanism in these surveys was lake water cooling. As the 
method proposed here is a first order estimate , it was 
decided for simplicity to ignore atmospheric cooling and 
utilize the relationship between isotheim and lake-plume 
temperature difference. 

Isotherm Distan ce and Temperature Difference 
(plume and ambient lake) 
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Although this relationship was not required 

to predict the isotherms, it is useful to obtain rapid 

estimates of maximum distances isotherms will occur away 

from the discharge point, without plotting actual isotherms. 

The equations are summarized in Table 4. It is observed 

that the equations were not as statistically significant 

as the other fits and the standard errors over standard 

deviation of distance were larger than the other fits. 

The best summer equation including a flow term is: 

Distance = 33.08 - 3,0 In (flow X 

temperature difference) (6) 

with F^ Qc- = 29.0: multiple correlation 

coefficient = 0.48: standard error = 2.18: 

standard deviation of distance = 2.48. 

The best winter equation including a flow term was: 

Distance = 26.88 - 2.30 In (flow X 

temperature difference) (7) 

with F, ,^ = 73.4: multiple correlation 
1,63 

coefficient = 0.7 3: standard error = 1.48: 

standard deviation = 2.16. 
The winter equation is the significantly better prediction 
equation (Figure 6) . 

Definition of Plume Centreline 

The curvature or radius of the plume was defined by 
fitting the measured radius of the plume centreline and the 
absolute value of the resolved component of the wind acting 
parallel to the shoreline. The locii for the centres for 
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the plume radii were determined by trial-and-error and 
appear in Figure 7. The generality of these locii in 
terms of other locations is not known particularly as 
there was an effect on the plume from the cooling water 
intake channel configuration. The best fit relationship 
for radii and wind component appears in Figure 8. 

R = 500 |wp| . (8) 

where R = plume centreline radius in ft. 

Wp= wind component along the shore in mph. 
It is apparent that the point scatter is large for equation 
(8) . The standard error of estimate can be improved by 
defining a unique locus of the centre of the radius for each 
experiment. However, this is of little value for a 
technique of representing the thermal plume. Other methods 
involving current measurements could not be used as the 
measurements are incomplete. 

Step Description of Method 

From the meteorological data, for the period under 
consideration, the wind component along the shoreline was 
computed and thence the radius of plume centreline from the 
equation (8). The plume centreline is drawn according to 
the definition of Figure 7. 

For the existing lake temperature, falling in 
'summer' or 'winter' regimes, the appropriate equations were 
selected from Tables 1 and 3. Plume temperatures and areas 
were computed at various distances from outfall and along the 
plume centreline for a particular discharge. For any 
distance, the area computed represented the plume contained 
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by the corresponding temperature isotherm. One leaf 
of the multiple-leaved rose was then fitted to the 
distance and plume area. A step by step computation 
appears in Appendix I. 

Three examples of the method are presented in 
Figures 9,10 and 11 for lake temperatures of 3 3 F, 46. 5° F 
and 71.7°F and wind vectors acting in three different 
directions. The actual field survey results appear as 
solid line isotherms whereas the predicted plumes are 
broken lines . The shaded area for the plume predictions 
represent one standard error of estimate for the plume 
area with isotherm enclosing the smaller area representing 
the mean value. Similarly, the temperature range is one 
standard error of estimate with the upper value representing 
the meam value. Errors of locating the plume centreline 
have not been included and the centreline represents the 
mean value location. 

DISCUSSION 

The isotherm predictions in Figures 9 to 11 do 
approximate the measured isotherms. The fact that the 
predictions are first order estimates with large standard 
errors is immediately apparent. Statistical fitting of 
data which knowingly do not encompass all the determining 
mechanisms, e.g. currents, thermal conductivity of lake, 
radiation, etc. has many shortcomings. In some cases, it 
would have been possible to estimate some of the mechanisms 
like currents and radiation from other surveys or serai- 
empirical equations. However, these estimations were 
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avoided as experience has indicated that currents very near to 
shore and near the surface tend to be very local in nature. 
Particular caution is also indicated in the Lakeview area 
where major shoreline modifications like the cooling water 
intake and wharf structure have taken place. 

The method is a first order estimate based on 
statistical fitting of limited data taken at one operating 
thermal generating station. It is suggested that the errors 
of estimate should be presented with the estimates to clearly 
identify the limitations of the methods. Transposing the 
predictions to another area has not been tested. 
CONCLUSIONS 

A method for predicting the isotherms of a surface 
cooling water discharge was developed statistically from 
temperature surveys of an existing thermal plume on Lake Ontario. 
However, the errors of estimating the plume were large indicating 
that a more comprehensive measuring program is needed to improve 
the predictions. The measuring program should be continuous 
over many hours and include measurements of all the important heat 
exchange mechanisms. The heat lost by mixing with colder lake 
water was significantly greater than the heat lost to the 
atmosphere (not considering evaporation implicitly). 
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TABLE 1 
Statistical Fitting 
Plume Area (A) and Distance from Outfall (X) 



Data 
Set 


Equation 


Std. 

Dev. 

A 


F 


F @ 

a=0.01 


Multiple 
Correlation 


Std. 
Error 


Summer V 


A* = -0.696 + 1.384X* 
A = 2.11 + 0.124X^ 
A = 4.38 + 0.00029e^ 


4.62 
4.62 
4.62 


F =117.1 

1 r 9 5 

F = 96.2 

1 f 9 S 

F = 28.5 

1^95 


6.91 
6.91 
6.91 


0.74 

0.71 
0.48 


3.11 
3.28 
4.07 


Winter^ 


A = -4.047 + 2.43X 5.80 
A = 0.91 + 0.239x^ 5.80 
A = 4.31 + 0.0033e^ 5.80 

1 1 


F =283.6 

1 f 6 3 

F =283.3 

1/63 

F =114.0 

1^63 


7.06 
7.06 
7.06 


0.90 
0.90 

0.80 


2.49 
2.49 

3.48 

1 



6 2 

Notes: # A (Area) in 10 ft. 



* X (Distance) in 10 ft. 

^ 'Summer' is defined when lake temperature >_45°F 

^ 'Winter' is defined when lake temperature <45''P 



TABLE 2 



Statistical Fitting 
Plume Area (A*) vs Temperature Differences AT ' (Plume -Lake) smd ^Tl,( Plume -Air) 



PL 



PA 



Data 
Set 


Equation 


Std. 

Dev. 

A 


W 


F @ 
a=0.01 


Multiple 
Correlation 


Std. 
Error 


• Summer' ^ 


A = 105.16 - 10.83 In (Q*AT ) 
+ 0.53 In (QATp^j ^^ 

QATpL : t = 10.65; t^ ^^=2.384 
Q%A = ^ = l-°3.' to 0^ =2.384 


4.676 


F =72.7 

2f 72 


4.91 


0.82 


2.73 


'Winter'^ 


A = 76.77 - 8.45 In (QAT ) 
+ 1.25 In (QATp^) 

QATp^: t = 8.98; t^ ^^ =2.395 

QATp^: t ^ 2.17; t^^^^ =2.395 


i 
5.194 


i 

F =42.6 

2, 57 

1 


5.00 


0.77 


3.35 



NOTES 



* A (Plume Area) x 10 ft^. 

# Q (Discharge from outfall) in cfs 



*# At (Temperature difference between plume and lake) F 

#* AT (Temperature difference between plume and air) F 

'Summer' is defined when lake temperature >45*^ 

2 O 

'Winter' is defined when lake temperature <45 F 



TABLE 3 



Statistical Fitting 
Plume Area (A) and (QATp ) 



Data 
Set 


-■ 

Equation 


Std. 
Dev. 

A 


F 


r ■ " ■'■ - 

F @ 

a= 0.01 


Multiple 
Correlation 


Std. 
Error 


' Summer ' * 


A* = 14.56 - 1.06 (AT**) 

A = 14.14 - 0.00055 
(Q*ATpj^) 

A = 23.67 - 8.79 In (AT„J 

PL 

A = 90.53 - 8.86 In (QATp^ 
A = 23.53 - 6.28 </ATpj^ | 
A = 23.25 - 0.14 >/QATp^ 


4.62 
4.62 

4.62 
4.62 
4.62 
4.62 


Fi, 95=112. 4 
Fi,95= 98.4 

Fi, 95 = 142. 5 
Fi, 95 = 135. 2 
Fi ,95 = 130.0 
Fi, 95=118. 9 


6.91 
6.91 

6.91 
6.91 
6.91 
6.91 


0.74 
0.71 

0.77 
0.77 
0.76 
0.75 


3.14 
3.26 

2.94 
2.99 

3.02 
3.10 


'Winter' ^ 


A = 15.91 - 0.89 (AT^^) 

PL 

A = 13.97 - 0.00034 (QAT^^) 
A = 23.70 - 7.81 In (ATp^) 
A = 74.00 - 6.88 In (QA"^) 
A =24.76-5.81 /AT^^ 

PL 

A =21.75-0.11 /QATpj^ 


5.80 
5.80 
5.80 
5.80 
5.80 
5.80 


Fi, 63 = 102. 9 
Fi ,63= 59.8 
Fi, 63=176. 8 
Fi, 63=126. 3 
Fi, 63 = 155. 4 
Fi,63= 94.3 


7.06 
7.06 
7.06 
7.06 
7.06 
7.06 


0.79 
0.70 
0.86 
0.82 
0.84 
0.77 


3.60 
4.19 
2.99 
3.37 
3.14 
3.70 



Notes : 



# 

* 



A (Area) in 10^ ft^ 

Q (Discharge) in cfs 

ATp (Temperature difference between plume and lake) in °F 

'Summer* is defined when lake temperature M5°F 
'Winter' is defined when lake temperature <45°F 



TABLE 4 

Statistical Fitting 

# 
Plrnne Distance (X) and Temperature Difference (AT )* 



Data 
Set 


Equation 


Std. 

Dev. 

X 


F 


F (a 

a= 0.01 


Multiple 
Correlation 


Std. 
Error 


' Summer ' ' 


X = 7.49 - 0.38 (ATpj^) 


2.48 


Fi, 9 5=29. 8 


6.91 


0.49 


2.17 




X = 7.39 - 0.00020 (Q#*ATpj^) 


2.48 


Fi , 9 5 = 28. 5 


6.91 


0.48 


2.19 




X = 10.25 - 2.9 In (ATp^^) 


2.48 


Fi, 95 = 27. 9 


6.91 


0.48 


2.19 




X = 33.08 - 3.0 In (QATp^^) 


2.48 


Fi, 95=29.0 


6.91 


0.48 


2.18 




X = 10.49 - 2.16 /ATpj^ 


2.48 


Fi, 95=29. 5 


6.91 


0.49 


2.18 




X = 10.48 - 0.05 /QAT^, 


2.48 


Fi, 95=29. 5 


6.91 


0.49 


2.18 


'Winter' ^ 


X = 7.73-0.32 (ATpj^) 


2.16 


Fi » 6 3= 87.8 


7.06 


0.76 


1.40 




X = 7.00 - 0.00012 (QATp^) 


2.16 


F 1 , 6 3- 51.3 


7.06 


0.67 


1.61 




X = 10.16 - 2.64 In (ATp^^) 


2.16 


Fa, 63= 98.4 


7.06 


0.7 8 


1.36 




X = 26.88 - 2.30 In (QATp^) 


2.16 


Fi,63= 73.38 


7.06 


0.73 


1.48 




X - 10.73 - 2.03 /ATpj^ 


2.16 


Fi, 63=106. 77 


7.06 


0.79 


1.32 




X = 9.61-0.04 *^QATpj^ 


2.16 


Fi,63= 67.78 


7.06 


0.72 


1.51 



Notes : 



#* Q 



X (Distance from outfall measured along plume centreline) in 10" ft 



AT 



PL 



(Temperature difference between plume and lake) in °F 



(Discharge from outfall) in cfs. 
'Summer' is defined when lake temperature >.4 5°F 
'Winter' is defined when lake temperature <45''F 



^^ 



r = X Cos n <E 




= 



^^ 



FIGURE I 
ISOTHERM SHAPE 
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APPENDIX I 



Sample Computations for Thermal Plume Predictions July 23, 1969 . 

Relevant Data 

Air temperature, T = 75°F 
Lake water temperature, T = 71.7°F 
Wind Speed, W = 10 mph @ 271° 
Outfall Discharge, Q = 1750 cfs 

Definition of Plume Centreline 

Since the shoreline is @ 231°, the angle between 
shoreline and prevailing wind is: 

= 271-231 = 40° 
. . Wind Component along shoreline 

Wp = W Cos e = 10 Cos 40 = 7.65 mph 
But R = 500 |Wp| ft 
. . R = 3830 ft. 

The arc of the plume centreline is drawn with a radius 
of 3830 ft. on locii of centres (see Figure 4). Since the 
arc does not lie wholly on the water, the plume centreline is 
smoothly joined to outfall, avoiding the shore. 

Definition of Plume Area 

Since the lake temperature is 71.7°F, the computations 
are done for 'summer' regime. The prediction equation is 
accordingly selected from Table 1. 

A = -0.696 + 1.384X ._ (1) 

fi y 

where A = area in 10 ft. 

X = distance from outfall, measured along 

-3 
the plume centreline m 10 ft. 

The following computations are for X = 4 (4000 ft.) 

(a) Mean Area 

A = -0.696 + 1.384(4) = 4.84 



(b) Maximum Area (One standard error) 

Standard error of equation (1) = 3.11 from Table 1. 

. . Maximum area A* = 4.84 + 3.11 = 7.95 

Determination of Temperature Difference i^'^pjj 

The prediction equation is selected from Table 3 as 
follows : 

A = 90.53 - 8.86 In (Q^Tp^^) (2) 

where ATp^ = Tp - T^ in °F 

Tp = plume temperature in °F 

(a) '^T (Mean Area) 

i: J_j 



Substituting the value A = 4.84, calculated above 
for X = 4 and Q = 1750 in equation (2) , ^Tp^ is computed. 

. . ATpj^ =9.05 

.'. T„ == T^ + At^^ = 71.7 + 9.05 = 80.75°F 
P L PL 

(b) At * (Maximum Area) 

XT Li 



Using A* - 7.95 and Q = 1750 in equation (2), 
At * is computed. 

At„ * = 6.37°F 

XT J-J 

.". T-* = T^ + At_ * = 71.7 + 6.37 = 78.07°F 
P L PIj 

Determination of Isothermal Contour 

Equation (3) describes a multi-leaved rose in polar 

co-ordinates (r,ct ) 

r = X Cos noi (3) 

where n governs the number of leaves and therefore the 

lateral spread. (See Figure 1) . 

Area enclosed by the curve may be shown as : 

A = 3x^ (4) 

4n 

For X = 4, substituting the values of A = 4.84 and 
A* = 7.9 5 in equation (3) values of n and n* are computed 

n = 2.6 

n*= 1.58 
Substituting these values of n and n* in equation (3) , the 
corresponding curves can be plotted. 









APPENDIX II 












Data Used 






Survey Date 


Discharge Outfall 
Flow Temp. 


Lake 
Temp. 


Air 

Temp. 






cfs 


°F 


% 


°F 


April 4, 


1969 


1970 


56.2 


40.0 


52.0 


May 22, 


1969 


1400 


62.0 


46.0 


50.0 


May 23, 


1969 


1400 


61.0 


46.0 


57.0 


June 19, 


1969 


1830 


58.6 


47.0 


60.0 


July 15, 


1969 


1730 


76.0 


61.7 


82.0 


July 16, 


1969 


1730 


78.0 


62.2 


83.0 


July 17, 


1969 


1620 


74.0 


61.7 


83.0 


July 22, 


1969 


1750 


84.0 


69.3 


75.0 


July 23, 


1969 


1750 


84.0 


71.7 


75.0 


July 24, 


1969 


1750 


86.0 


71.5 


76.0 


Sept. 9, 


1969 


1970 


72.0 


59.0 


66.0 


Sept. 10 


, 1969 


1750 


70.0 


58.0 


69.0 


Sept. 11 


, 1969 


1750 


74.0 


59.0 


63.0 


Sept. 30 


, 1969 


1830 


73.0 


58.0 


66.0 


Oct. 1, 


1969 


1830 


73.0 


5 8.0 


60.0 


Oct. 2, 


1969 


1980 


75.0 


59.5 


61.0 


Oct. 21, 


1969 


1740 


62.0 


43.0 


48.0 


Oct. 22, 


1969 


1760 


58.0 


43.0 


31.0 


Oct. 23, 


1969 


1720 


56.0 


42.0 


33.0 


Nov. 12, 


1969 


2130 


64.0 


47.0 


47.5 


NOV. 13, 


1969 


2140 


60.0 


46.5 


42.5 


Nov. 25, 


1969 


2040 


57.0 


41.0 


45.0 



APPENDIX II (Cont'd) 
Data Used 



Survey Date 


Discharge 
Flow 


Outfall 
Temp. 


Lake 
Temp. 


Air 
Temp. 






cf s 


°F 


^P 


°F 


Nov . 26, 


1969 


2030 


56.0 


41.0 


40.0 


Nov. 27, 


1969 


2010 


56.0 


41.0 


32.5 


Jan. 28, 


1970 


2520 


56.0 


35.0 


36.0 


Feb. 3, 


1970 


2190 


51.7 


3 4.0 


8.0 


Feb. 4, 


1970 


2150 


51.4 


3 3.0 


13.0 


Feb. 20, 


1970 


2205 


61.0 


35.0 


22.0 


May 20, 


1970 


1924 


63.0 


45.0 


52.0 


May 21, 


1970 


192 


64.0 


46.0 


55.0 


June 10 , 


1970 


18 6 


66.0 


61.0 


71.0 


June 11 , 


1970 


1850 


66.0 


54.0 


71.0 


June 11, 


1970 


1740 


65.0 


54.0 


84.0 


Sept. 1, 


1970 


2270 


65.0 


56.0 


65.0 


Sept. 2, 


1970 


2440 


68.0 


59.0 


60.0 


Sept. 2, 


1970 


2380 


67.0 


55.0 


65.0 


Oct. 20, 


1970 


2080 


70.0 


54.0 


53.0 


Oct. 20, 


1970 


1980 


72.0 


54.0 


54.0 


Oct. 22, 


1970 


1970 


72.0 


55.0 


55.0 


Oct. 22, 


1970 


1630 


72.0 


55.0 


56.0 



(16019) 
MOE/LAK 0/EST/APEN 



MOE/LAK 0/EST/APEN 
Palmer, M D 
Estimate of the 
physical extent of 
rtemu' piuinci "^ x 
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